,  yftD-A13-4  427 


UNCLASSIFIED 


HOMOGENIZATION  FOR  A  VOLTERRA  EQUATIONS)  WISCONSIN' 
UNI V-HADISON  MATHEMATICS  RESEARCH  CENTER 
H  ATTOUCH  ET  AL,  SEP  82  MRC-TSR-2566  DAAG29-80-C-0041 

F/G  12/1 


MRC  Technical  Summary  Report  #  2566 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53705 


September  1983 


(Received  March  2,  1983) 


Approved  for  public  release 
Distribution  unlimited 

NIC  FILE  COPY 

Sponsored  by 

U.  S.  Army  Research  Office 

P.  O.  Box  12211 
Research  Triangle  Park. 

North  Carolina  27709 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 


HOMOGENIZATION  FOR  A  VOLTERRA  EQUATION 

*  ** 

Hedy  Attouch  and  Alain  Damlamian 

Technical  Summary  Report  #2566 

September  1983 

ABSTRACT 

A  model  is  given  for  the  non-linear  heat  equation  in  a  heterogeneous 
medium  with  memory.  Its  homogenization  is  carried  out  in  two  particular  cases 
(including  the  linear  one). 
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SIGNIFICANCE  AND  EXPLANATION 


Nonlinear  heat  flow  in  a  heterogeneous  material  is  considered.  In  this 
model,  the  Internal  energy  and  heat  flux  depend  upon  the  history  of  the 
temperature  and  the  gradient  of  the  temperature  respectively. 

The  heat  conservation  law  leads  to  a  Nonlinear  Vo 1 terra  integro- 
differential  equation  with  appropriate  boundary  conditions.  This  problem  is 
solved  under  physically  reasonable  assumptions  and  its  homogenization  is 
investigated:  introducing  a  small  parameter  e  measuring  the  "tightness"  of 
the  heterogeneity  of  the  medium  (typically  we  assume  e-periodicity  for  the 
physical  parameters ) ,  the  stability  of  the  model  is  studied  (as  e  goes  to 
zero)  and  the  homogenized  (ideal)  limit  medium  is  characterized  in  some  cases, 
including  the  linear  one. 


The  responsibility  for  the  wording  and  views  expressed  in  the  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


HOMOGENIZATION  FOR  A  VOLTERRA  EQUATION 


*  ** 

Hedy  Attouch  and  Alain  Daralamian 

I.  Introduction 

In  a  heterogeneous  medium  with  memory,  a  model  for  the  heat  equation  (see 
Nohel  [1])  is 

(1.1)  |^  +  divx  Q  -  h1  , 

where  h^  is  a  given  diffused  source  term,  €  is  the  internal  energy,  Q  is 
the  heat  flux.  The  latter  are  assumed  to  be  functionals  of  the  temperature 
distribution  u  with  "memory": 

(1.2)  S(t,x)  «  bQ(x)u(t,x)  +  8(x,  t-s)u(s,x)ds 

Q(t,x)  -  -cQ  (x)a(x,Vu(t,x) )  + 

(1.3) 

+  Y(x,t-s)a(x,Vu(s,x) )ds 

Equation  (1.1)  is  considered  on  the  product  0  x  (-»,T),  where  0  is  a 

J  M  J|| 

bounded  regular  domain  in  R  (or  Rr)j  the  function  a  :  Si  x  r  -*■  R  , 

(x,r)  ♦  o(x,r)  represents  a  nonlinear  flux  law.  Its  dependence  upon  x 
specifies  the  heterogeneity  of  the  medium.  Similarly  bg(x),  6(x,t),  cQ(x), 
Y(x, t)  characterize  the  spatial  heterogeneity  of  the  other  thermodynamical 
parameters. 

TO  equation  (1.1)  are  added  boundary  and  initial  conditions  which  will  be 
specified  later. 

The  questions  considered  here  are: 

'  * 
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-  under  what  suitable  set  of  hypotheses  is  equation  (1.1)  well-posed 


(existence  and  uniqueness); 

-  under  what  further  conditions  can  one  treat  the  corresponding  homogenization 

problem;  in  other  words,  if  all  the  parameters  involved  (c,bg ,B,Cg ,y)  depend 

upon  another  variable  e  measuring  the  "tightness"  of  the  heterogeneity  of 

the  medium  (typically  ^(x)  *  bg(“)  where  b(j(y)  is  periodic),  can  one  find 

e 

a  limit  problem  whose  solution  would  be  the  limit  of  the  solutions  u  ,  and 
whose  structure  would  be  similar  to  (1.1),  (1.2),  (1.3)7 

We  will  give  positive  answers  to  both  questions  in  particular  cases 
only:  in  one  case  the  nice  method  of  Crandall-Nohel  [1]  applies  and 
homogenization  follows  (§3,  the  splitting  case);  in  §4,  we  deal  with  the 
general  situation  for  which  existence  and  uniqueness  is  proved  via  local 
monotonicity  and  global  estimates;  in  $5  we  treat  the  homogenization  for  the 
linear  case  via  the  Laplace  transform. 

A  first  approach  to  this  type  of  problem  appeared  in  Raynal  [1], 

II.  Reformulation  of  the  problem  as  a  Volterra  equation. 

Let  *  denote  the  usual  convolution  with  respect  to  t  on  [0,+«®[. 

As  "initial"  condition,  we  assume  the  history  of  the  medium  to  be  known 
for  t  negative.  One  can  then  rewrite  (1.1),  (1.2),  (1.3)  as: 

(2.1)  —  [bQu  +  f3*u]  +  divx(-c0o  +  y*o )  -  h 

where  a  stands  for  o(«,Vu(*,»)  and  the  right  hand  side  h  includes  the 
history  of  the  system  up  to  time  zero. 

It  is  customary  to  define 

(2.2)  c(x,t)  =  cQ(x)  -  /q  Y(x, s)ds 

and  to  assume  that  c(x,t)  and  bQ(x)  are  strictly-positive  valued  (a 
physical  condition).  With  these  notations  (2.1)  can  be  written  as: 


(2.3) 


bgu'  -  divx{ (c#o) ' }  ■  h  -  8gU  -  f5’*u 


t 


where  indicates  a  time  derivative. 

The  initial  condition  becomes  a  Cauchy  data  at  t  *  0: 

(2.4)  u(0)  -  uQ  in  L2(fl)» 

the  boundary  condition  is  taken  to  be  compatible  with  the  operator 
-div^(o( Vu) ) ,  for  example 

(2.5)  u  -  0  on  3fl  x  (0,T)  . 

Problem  (2.3),  (2.4),  (2.5),  upon  integration  with  respect  to  time 
appears  as  a  Volterra  type  integral  equation  (cf.  3.1). 

We  now  make  precise  the  type  of  function  o  appearing  here:  let 
j  :  (x,r)  j(x,r)  e  *+,  be  of  Caratheodory  type,  convex  and 

equicoercive  in  r.  Assume  further  that  j(x,0)  =  0;  let  a  be  its  sub¬ 
differential  with  respect  to  r  and  put: 


(2.6) 


4(u) 


/ n  j(x,Vu(x))dx  for  u  e  w^lfl) 
+*>  otherwise  i 


then  it  is  shown  in  Attouch-Damlamian  [1]  that  $  is  lower  semi-continuous 

2 

(l.s.c.)  convex  on  L  (fl)>  moreover,  denoting  by  3$  the  subdifferential  of 

2 

4  on  L  ( ft )  one  has : 


v  e  3$(u) 


u  e  Wg#1(n)  o  l2(Q) 


\  v  ■  -div  h  for  some  h  such  that 
h(x)  e  o(x,7u(x))  a.e. 

This  is  the  sense  in  which  a  is  used  in  (2.3)  and  it  also  gives  a  meaning  to 
the  Dirichlet  condition  (2.5)  (when  j  is  even  with  respect  to  r,  (2.7)  is 


actually  an  equivalence) 


is  independent  of  x.  Equivalently,  by  an 


III.  The  Splitting  Case 

c 

By  this#  we  mean  that  — 

co 

obvious  change  of  notation  (for  a),  c  and  Y  can  be  taken  independent  of 
x  (cfl  is  taken  to  be  1 ) . 

A)  Existence  and  Uniqueness 

Integration  of  (2.3)  with  respect  to  time  from  0  to  t  yieldsi 

(3.1)  bQu  -  divx(c*o)  *  b0u0  -  0*u  +  H 
where  R  is  the  integral  of  h. 

—1  »  + 

(3.2)  Proposition:  Assume  that  bQ  and  bQ  are  in  L  (0)  ,  that  B 

is  in  L  (ft;  BV(0,T)),  and  that  Y  is  in  BV(0,T).  Then  equation  (3.1) 

2  1  11 

has  a  unique  solution  u  in  C( [0,T] ;  L  (ft))f\  L  (0,T,WQ'  (ft)), 
du  2 

Furthermore,  ~  is  in  L  (ft  *  (0,T)). 
at 

Proof :  we  follow  here  the  ideas  of  Crandall  -  Nohel  [1].  Let  e  be  the 
resolvant  kernel  of  Y,  i.e.  the  solution  of 

(3.3)  e  -  Y  -  Y*e  -  0  . 

Making  use  of  standard  results  for  convolution  equations,  one  obtains 
that  e  belongs  to  BV(0,T)  as  soon  as  Y  does  so.  Using  (3.3),  (2.3) 
becomes 


(  bQ  —  -  divx  o(x,Vu)  -  G(u) 

^  u(0)  -  uQ  ,  where 

G(u)  «  h  +  h*e  -  (B  +  b  e. )u  +  bue 

0  0  0  0  0 


(3.5) 


-  u  *  (0Qe  +  bQe'  +  B*  +  e*B';) 

Above  Bq  ”  B(0),  e^  ■  e(0)  (■  Y(0))  and  B',  e';  are  the  measure 
derivatives  of  B  and  e. 

It  is  easy  to  check  that  G  is  Lipschits  continuous  from 


1  2 

L  (0,t;L  (ft))  into  Itself  for  each  t  >  0.  In  order  to  apply  a  fixed 


point  theory  (as  in  Crandall  -  Nohel  [1])  we  first  consider  the  following 

2 

problem  for  w  in  C{  [0,T] ;  L  (Jl)): 

{b  ~  -  div  o(x,Vu)  ”  w 
0  d* 

u(0)  -  Uq  . 

Now,  the  operator  u  ♦  -  ““  div  a(x,Vu)  is  the  subdifferential  of 

*  D0  x  . 

2  2 

4  (see  (2.6))  on  L  (0)  provided  the  norm  on  L  (ft)  is  chosen  with  the 

00  + 

weight  function  bg ( x)  (since  bg  and  bg  are  in  L  (ft)  ,  this  is  an 

2 

equivalent  Hilbert  norm  on  L  (ft)).  Therefore,  (3.6)  can  be  solved  with 
classical  estimates  which  allow  to  apply  the  Lipschitz  fixed  point 
theorem  to  solve  (3.4),  (3.5).  ■ 

B]  Homogenization 

e  e  e  e  c 

We  now  assume  that  bg ,  cQ ,  Y  ,  B  and  a  depend  upon  an  extra 

parameter  e  measuring  the  size  of  the  heterogeneity  of  the  medium.  A 

0  x  E 

typical  example  is  the  periodic  case  where  bg(u)  *  bQ (— ) ,  cQ(x)  - 
Cg(~),  etc.  ...  where  bQ(y),  Cg(y),...  are  Y-periodic  (Y  is  an  N- 

dimensional  parallelepipedon) .  We  make  the  following  hypotheses: 

E  E  -1  oo  + 

bg  and  (bg)  are  bounded  in  L  (ft)  , 

$E  is  bounded  in  L  (ft»  BV(0,T))  and  yE  i*  bounded  in  BV(0,T). 

Applying  proposition  (3.2)  one  gets: 


Proposition:  Under  hypothesis  (3.7),  there  exists  a  unique  solution 

£  £  2 
u  for  problem  (3.1)£  and  u  is  bounded  in  C([0,T]»  L  (ft))  by  a 


constant  involving  only  |h|  2  ,  1(5-1 


.dlnv  '  "O'  ®  '  ^0^  °»  ' 

L  (Q)  w  L  (ft)  L  (0) 


I  *»0>— 1 1 


L  (0) 


,  I6e  I 


L  (0;  B V { 0 , T ) ) 


*"'1  |T  ^BV( 0,T)  * 


e 

In  order  to  study  the  convergence  of  u  when  e  goes  to  zero  we  make 

the  following  extra  hypotheses  (3.9)  -  (3.13): 
e 

(3.9)  j  (x,r)  is  coercive  in  r  uniformly  with 

respect  to  x  and  £;  more  precisely: 

2N 

a  >  0,  3  K2  e  R  and  p  >  such  that 

2N 

.£,  >  .  v  I  .N+2 

j  (x,r)  >  K1  1 r |  -  *2  , 

£  (1)2 

(3.10)  $  (defined  as  in  (2.6))  converge  in  the  sense  of  Mosco  on  L  (fl) 

to  a  limit  denoted  $,  which  is  then  known  to  be  of  the  same  integral 
form,  associated  to  a  convex  function  j  (cf.  Attouch  [1]).  Finally 

t(uQ)  is  assumed  to  be  finite. 

£  oo 

(3.11)  P  converges  in  the  weak-star  topology  of  L  (fl»  BV(0,T))  to  a  limit 

e 

P  (consequently  PQ  converges  to  PQ  in  the  weak  star  topology  of 

■  c 

L  (£));  in  the  periodic  case,  P  is  just  the  average  of  P  over  a 

period) . 
c 

(3.12)  bQ  converges  to  some  b0  in  the  weak  star-topology  of  L  (Q)  and 

eE  converges  to  some  e  in  the  weak-star  topology  of  BV(0,T). 

Then,  e  is  the  resolvant  of  some  Y  in  BV(0,T)  (but  Y  has  no 

relationship  to  the  weak-star  limit  of  YE  in  BV(0,T)). 

£ 

Consequently,  the  mapping  6  (the  analogue  of  6  in  (3.5))  is  bounded 
e  £  2  2 

so  that  G  (u  )  is  bounded  in  L  (0,T:  L  ((J))  and  therefore,  via  the 
properties  of  the  solutions  of  the  problems  (3.6)e,  one  can  conclude  (making 
also  use  of  (3.10))  that 


(1) 

For  the  definition  and  properties  of  this  convergence,  see  Mosco  [1]  or 
Attouch  [2]. 


is  bounded  in 


(3.13)  u  is  bounded  in  L  ( 0 , T ;  Wn  ( ft ) ) ,  - — 

o  at 


2  2 

L  (0 , T>  L  (.SI)),  hence,  as  a  consequence  of  Aubin's  lemma  (see  Aubin 


£  2 
[1]),  u  is  compact  in  C([0,T}>  L  (0)). 


We  will  show  now  that  u  has  only  one  possible  limit  value  u  when  £ 


goes  to  zero,  which  implies  that  u  converges  to  u.  Let  therefore  £ 


go  to  zero  so  that  u  n  converges  to  some  u  in  C(  [ 0 , T ] f  L^{£2)): 


e  e 

(3.14)  Proposition:  Under  the  above  hypotheses  ((3.13)),  G  n(u  n) 


»nd 


converges  weakly  in  L  (£2  x  (0,T))  to  G(u)  (as  given  in  (3.5)),  and 


£  .  n  2  j 

n  du  au 

b„  — — —  converges  weakly  in  L  (ft  x  (0,T))  to  b„  — . 
0  dt  0  dt 


Proof:  We  write  £  instead  of  £  for  simplicity,  bf  is 

-  n  0  dt 


bounded  in  L  (Q)(Q  =  £2  x  (0,T))  and  for  <p(x,t)  in  V  (Q)  one  has: 


/  be  =  ./  be  if  b  if  u 

JQ  0  dt  JQ  0  at  JQ  0  3t 


e.  £. 


which  proves  the  second  claim.  For  G  (u  ),  each  term  can  be  treated 


independently;  let  us  show  convergence  for  the  worst  case; 


.£'  E^  £ 


8  *u  *e  .  First 


wt(x,t)  =  /q  eE(t-s)  uS(x,s-a)d8t'(x,o)ds 


so  that 


|wE(x,*)|  „  <  I e E |  IuE(x,*)|  ^  I 0E(X, * ) lgv(o  T)  ' 

L  ( 0,T)  L  ( 0,T)  L  ( 0,T)  ' 


<  |eE| 


iuE| 


8 


L  (Uf  L  ( 0,T ) ) 


L  ( 0,T )  L  (£2;  L  (0,T))  L  (£2;  BV(0,T)) 


£  2  00  £ 
Consequently,  w  is  bounded  in  L  (£2;  L  (0,T))  because  u  being 


2  1  2 

bounded  in  L  (£2;  H  (0,T))  is  bounded  in  L  (£2;  C([0,T])).  By  a 


similar  argument,  one  checks  that  for  almost  every  x  in  £2 


e 

u  (x,t)  >  u(x#t)  in  C( [0  ,T] )  which  will  be  enough  to  show  the 

e  rw 

convergence  of  w  to  the  proper  w  in  P*  (Q)  as  follows:  for  * 
in  P(Q) 

<wE^>>  -  ”/jj/J/  <?' (x,t+s+a)eC(t)ue(x,s)Be(x,o)  dt  ds  do  dx 
*  //  dtds  e6(t)//jj  do  dx  (x,t+s+o)ue(x,s)Be(x,o)  . 


For  every  (t,s)  the  last  integral  converges  to 

/  do  /q  V>'  (x,t+s+o)u(x,s>f5(x,o)dx 
e  a 

because  0  converges  to  3  *— weakly  in  L  (Q).  Furthermore,  Lebe  ■) 1  s 

dominated  convergence  theorem  applies  to  the  (t,s)  integral  since 

integrand  is  actually  bounded  by  a  constant,  namely 

l*'L  *  sup|ue|  w  •  sup | PE|  , 

L  (0,T»  L  (ft))  L  (Q) 

c 

this  last  factor  being  bounded  above  by  sup | 3  |  which  is  finite 


L  (ft,BV(0,T) ) 


by  hypothesis  (3.7).  ■ 


e 

Making  now  use  of  the  convergence  in  the  sense  of  Mosco  of  <p  to  <p 
which  implies  a  demi-closedness  property  (cf.  Attouch  [2])  one  passes  to  the 
limit  in 


-div  oe(x,Vu6)  ■  G£(ue) 

.  e  du 
“  bo  dt 

to  conclude  that  u 

is  a  solution  of 

(3.15) 

C  -div  o(x,Vu)  ■  G(u)  - 

bo  dt 

l^u(0)  -  uQ  . 

Because  G  is  as  in  (3.5)  one  concludes  to  the  uniqueness  of  the 
solution  for  (3.15),  hence  the  conclusion: 

(3.16)  Theorem:  Under  the  hypotheses  (3.7),  (3.9),  (3.10),  (3.11),  (3.12) 

e  e  E 

(and  with  the  notations  therein),  the  solution  u  of  problem  (3.4)  ,  (3.5) 


-8- 


converges  uniformly  in  C{[0,T]»  L  (A))  to  the  solution  u  of  the  analoguous 

problem  whose  thermodynamical  parameters  are  obtained  as  follows : 

€  £ 

0  and  bQ  are  the  weak-star  limits  of  0  and  bQ 

respectively 7 

e 

0  is  the  elliptic  homogenization  of  a  (equivalent  to  the 
Mosco  convergence  of  $E  to  <p  c f.  (3. 10)  >7 

e 

y  is  the  resolvent  of  the  weak-limit  of  the  resolvant  of  y 
( these  two  operations  do  not  commute  I ) .  ■ 


IV.  Existence  and  uniqueness  in  the  general  case 
He  start  with  equation  (2.3)  again 

(4.1)  bQu'  ”  div((c*o)')  »  g(u),  u(0)  ■  u^  , 

where 


(4.2)  g(u)  -  h  -  0Qu  -  0  '*u 

1  2 

is  Lipschitz  continuous  from  L  (0,t;  L  (fl))  into  itself  for  every  positive 
t.  Therefore,  one  can  solve  (4.1)  as  a  Lipschitz  perturbation  problem  (in  a 
fashion  similar  to  that  in  XII  A) )  provided  one  can  solve 


(4.3) 


bQ  ~  -  div( (c*a ) ' ) 


F 


u(0  )  -  uQ 

for  given  F,  via  a  monotonicity  argument. 

Here,  one  should  notice  that  the  method  of  Crandall  -  Nohel  does  not 
aPPly  because  c  depends  upon  x  so  that  div  and  convolution  with  c  Jo 
not  commute. 

He  make  the  following  assumptions  where  a  and  k  are  given  positive 
numbers : 


i) 


1  00 
0n,  t—  are  bounded  by  k  in  L  (ft)  and 


(4.4) 


1 


ii)  8'  <  0  and  c'  ■  -Y  <  0  a.e.  in  Q;  8(T,x)  and  c(T,x)  are 
bounded  below  away  from  zero  by  oj  t*+  c(t,x)  is  continuous 
at  t  ■  0  with  values  in  L  (ft). 

iii)  8"  and  c"(“  -Y ' )  are  nonnegative  measures  for  a.e.  x  in 

0. 

iv)  (2,6)  holds  with  the  following  inequalities 

V  r,s  in 
2 

a|r-8|  <  (o(x,r)  -  a(x, s),  r-s) 

and  |<J(x,r)  -  9(x,s)j  <  k|r-s|  . 

We  start  by  choosing  t  small  enough  as  follows: 

(4.5)  Proposition:  Let  Au  *  -div((c*o)')  *  -div(co0)  -  div(c'*o).  For  T 

2  1 

small  enough,  A  is  maximal  monotone  from  L  (0,T»  H^((2))  «  V  into 
L2(0,T);  H-1)  “  V'. 

Proof:  We  estimate 
(Au-Av,u-v)v,  v  > 

>  Jq  cQ{o(x,Vu(x)  -  o(x,Vv(x)),  Vu(x)  -  7v(x))dxdt 

+  jg  (c'*(o(x,Vu(x) )  -  o(x,7v(x)),  Vu(x)  -7v(x))dxdt 

>  a  /  c„ | 7u-7v| 2dxdt 

1 Q  0 

-  /g  )c-l*|o(Vu)  -  o(7v)|  |7u-7v|dxdt 

>  a  / g  cQ|Vu-Vv|  dxdt 

-  Jn  dx|c'(x)|  .  |a(7u)-o(7v) |  |7u-7v| 

“  L  (0,T)  L  (0,T)  L  (0,T) 

>  a  L  c.|7u-7v|2dxdt  -  k|c*|  m  ,  J  |7u-7v|2dxdt  . 

U  L  (£2>  L  (0,T) )  U 

Note  now  that  |c'(x)|  .  =  cn(x)  -  c(t,x)  so  by  (4.4)  ii) 

L1 ( 0,T )  0 

|c'|  ^  ,  is  arbitrarily  small  for  T  arbitrarily  close  to  zero. 

L~(fl|  L1(0,T)) 

2 

For  such  a  T,  (Au-Av,u-v)  >  9|7u-7v|  for  some  positive  number  9. 

L  (Q) 
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Hence  A  is  monotone.  Being  aiso  everywhere  defined  and  continuous  on  V, 
it  is  maximal.  ■ 


Now,  for  such  a  small  T,  brt  ~  +  A  is  one  to  one  and  onto  from  V  to 

o  at 

V'  (because  of  a  standard  non-linear  argument  of  coerciveness,  cf.  Brezls 
[1]).  This  proves  local  existence  and  uniqueness  of  the  solution  for  problem 
(4.3),  hence  for  (4.1),  (4.2).  In  order  to  prove  global  existence  we  now  get 
two  a  priori  estimates. 

(4.6)  Proposition:  Under  the  above  hypotheses,  there  is  a  constant  C., 
(depending  upon  a  and  k)  such  that  if  u  is  a  solution  of  (4.1),  (4.2), 
the  following  holds: 

I u |  ^  .  <  C  and  |j(Vu)|  ^  ,  |j  (o(Vu))|  .  <  C 

L  (0,T»L  (fi))  L  (Q)  L  (Q> 


( consequently  | Vu |‘ 


2  <  kC1 ,  |o(Vu) |  2 


I*  (Q) 


L  (Q) 


Proof :  Multiply  (4.1),  (4.2)  by  u(t)  and  integrate  by  parts  to  get 

Ib02u(t)|2  +  l^»(t)|22 

2  dt  °  L2(Q)  °  L2(0) 

(4.7)  +  /^(P'*u) (t)u(t)dx  +  f^(cQo+c'*a)(t)  •  Vu(t)dx  - 

=*  / q  h(t)u(t)dx  . 


Integrating  on  (0,t)  yields 


J  |bp2u(t)|22  +  /q  /n  BQu2dsdx  +  /*  fa  cQo(  7u(  s)  )Vu(  s)dsdx 


L‘(0) 


(4.8) 


■  / p  h»u  dsdx  +  /J  /n  (-f5,#u)  (s)u(s)dsdx  + 

+  /«  /0  (-c'*o)(s)Vu(s)dsdx  +  ^  lb^ u_| 2_ 
on  2  !,  (SI) 


By  (4.4)  ii),  one  has 
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(4.9  )< 


/0(-S'*u)(s)u(s)ds  <  J*  J*  -3* (s-T) (-|u(s) |  +  — I u( T ) |  )dsdT 

<  J|u(s)|2(30-6(8))d8  +  J  /q  da  /*  -$'(8-T)|u(T)|2dT 

i  <  £  /J(60-P(8))|u(8)|2d8  +  “  (/J  -3'(T)dT)(/J|u(T)|2dT) 


2  '•'O 


<  (BQ-3(t))  /q  u  (s)ds  . 


Hence 


(4.10)  J*  /jj  SqU2  +  (S’*u)u  dxds  >  /jj  B(x,t)  /q|u(x,s)  1 2dsdx  . 

A  similar  computation,  making  use  of  Young's  inequality  corresponding 
* 

to  j  and  j  ,  gives: 


(4.11) 


!l  /jj(co0(Vu)  +  c'*o(Vu))  •  7  u  dxds  > 

>  /n  c(x,t)  /q  j(x,7u(x,s))  +  j  (x,0(x,Vu(x,s) )dsdx  . 


Now  (4.8),  (4.10)  and  (4.11)  combined  give 

,  |b^2u(t)|2  +  L  8(t,x)  j!j  u2(s,x)dsdx  + 

L  (fl) 


(4.12) 


+  /n  c(t,x)  /g(j(s,Vu)  +  j  (x,0(Vu)))dsdx  < 

<  J$ /„  h  «  ^  I^OV,.  . 


v  '°'0 -  2  '  °  0  l2<o> 

A  standard  application  of  Gronwall's  inequality  finally  yields  the  desired 
result. 


(4.13)  Proposition:  Assume  the  above  hypotheses  and  that  4(Uq)  is 
finite.  Then,  there  is  a  constant  C2  (depending  upon  a,  k  and  $ ( uQ ) ) 
such  that  whenever  u  is  a  solution  of  (4.1),  (4.2)  then: 

lf?l  2  <  C2  '  !*(u)|  -  <  C2  * 

Z  L  (Q)  L  (0,T) 


fill 

Proof :  Multiply  (4.1),  (4.2)  by  —  to  get 


6 


/J  4  bo,d?|2dxd#  +  /S  Veou  +  e'*u>  (8)dxd8  + 

(4.14) 

+  fo  ^fl(cO°  +  c'*0)  dt  {7u<*>)dxds  *  Jo  Jn  h  ^(«)dxds  . 

In  (4.14)  we  Integrate  by  parts  the  third  term  as  follows: 

Jj  (e,*u)u*  -  u(t>  Jj  u(s)8' (t-s)ds  -  3 * ( 0 )  /J  u2(s)ds 

-  /q  u(s)  /8  u(o)$"(s-o)dods  . 

Making  use  of  8"  >  0,  one  can  evaluate  the  last  term  in  a  way  exactly 
similar  to  (4.9)  above  to  get 

(4.15)  /J  (8'*u)  >  u(t)  /J  u(s)8* ( t-s)ds  -  3* (0 )  /J  u2(s)ds  . 

Again,  similarly, 

(4.16)  /b  (c ,#cr)  “  (Vu)ds  >  7u(t)  •  (c**o)(t)  -  c'(0)  /b  o(s)7u(s)ds  . 
Using  (4.15),  (4.16)  and  the  following  consequence  of  the  definition  of  sub* 
differential  a  -  3j: 

(4.17)  0  •  (7u)  -  j(7u>  , 

(4.14)  yields 

bn  lSl2dxd®  +  2  l^®(t>|22  + 

0  dt  20  I#2(fl) 

/  cQj(Vu(t))dx  <  Iq  Ia  *»  dxd®  + 

2  ^  V,,V*‘  * 

u  \  U  I 

/n  u(x,t)  /q  u(x,s) (-5* (x,t-s) )dxds 


f*  I 
Jo  Jn 

+ 

(4.18)  + 


+  7u(x,t)  Jb  o(x,7u(x,s) ) t-c’ (x,t-s) )dsdx 


In  the  right  hand  side  of  (4.18)  one  can  use  the  following  bounds: 
i)  Jn  /q  u( t)u(s) (-8 ’ (t-s) )dxds  4 

(4.19) 


<  lul 


IS'  I 


L  (0,T»L  (fl)) 


L1(0,T;L*(n)) 


13 


ii)  By  Young's  inequality 


/jj  /q  o(s)Vu(t)  •  (-o'  (t-s)  )dsdt 

<  /n  j(Vu(t))(c.-c(t))  +  |c'|  m  fj  (o(Vu))|  . 

It  (Q)  L  (Q) 

Mow,  confronting  (4.18)  and  (4.19)  with  hypothesis  (4.4)  ii)  and  proposition 
(4.6)  one  can  conclude.  ■ 


(4.20)  Theorem:  Under  hypotheses  (4.4),  problem  (4.1),  (4.2)  has  a  unique 
solution  on  [0,T]. 

Proof:  From  proposition  (4.5),  there  is  existence  and  uniqueness  on  some 
interval  [0,t]  t  >  0.  But  the  very  same  proposition  gives  existence  and 
uniqueness  locally  in  time  (starting  from  t  >  0,  the  problem  is  changed  only 
insofar  as  G  in  (4.2)  is  modified  to  incorporate  the  history  up  to  t;  this 
in  no  way  changes  the  conclusion  because  of  the  a  priori  estimates  are  global 
in  time).  Combining  local  existence  and  the  a  priori  estimates  (4.6)  and 
(4.13)  gives  the  result  in  a  standard  way.  ■ 


V.  Homogenization  (for  the  linear  case) 

C  C  C  £ 

In  this  paragraph,  bQ,  3  ,  c  and  y  will  depend  upon  the  parameter 

£ 

e  which  will  tend  to  zero:  similarly,  a  will  depend  upon  e  but  will  be 

assumed  to  be  linear  with  respect  to  Vu,  hence  the  notation 

£  £  £ 
o  (x,r)  *  A  (x)r  ,  where  A  (x) 

is  a  measurable  function  from  0  to  a  fixed  (independent  of  e)  set  of 
symmetric  uniformly  positive  definite  matrices. 


We  shall  assume  hypotheses  (4.4)  to  be  satisfied  uniformly  with  respect 


to  e,  so  that  (4.6),  (4.13)  and  (4.20)  hold  uniformly  in  e. 

Consequently,  the  solutions  u6  of  equations  (4. 1)  ,  (4.2)e»  belong  to 

2 

a  compact  set  of  C([0,T]r  L  (0))  and  a  bounded  set  of 


L  ( 0,T;  HQ(fl))  O  W  '  (0,T»  L  (ft))  . 

The  question  of  homogenization  for  (4.1)  ,  (4.2)  is:  what  can  be  said 

e  e 

e 

of  the  limit  points  of  u  as  e  goes  to  zero? 

In  order  to  simplify  the  notations  we  shall  assume  that  e  belongs  to  a 

sequence  (still  denoted  c)  such  that  the  following  holds: 

E  2 

(5.1)  u  converges  to  some  u  in  C((0,T]»  L  (0)),  in  the  weak-star 

00  t 

topology  of  L  (0,T;  Hg(fl))  and  the  weak  topology  of 
W1,2(0,T;  L2(fl))i 

£  oo  0 

(5.2)  bQ  converges  to  bg  in  the  weak-star  topology  of  L  (fl);  B 
converges  to  B  in 

o(W  1#1(0,T»  l“(Q)),  W_1'"(0,T»  L1  (Jl) ) )  . 


Following  integration  in  t,  the  Volterra  equation  can  be  written  as 

(5.3)  -  div  W6  »  F6  , 


where 


(5.4) 


W£(x,t)  ■  AE(x)c£(xrt-s)Vue(x,s)ds 


£  £  £  f  0  0 

F  (x,t)  *  -bQ(x)u  (x,t)  -  / Q  B  (x,t-s)u  (x,s)ds  + 

(5.5) 

+  / g  h(x,s)ds  +  bg(x)uQ(x) 

_£  -1 
Clearly  r  converges  to  r  in  C((0,T],  H  (0))  and  weakly  in 

w1,2(0,Tj  L2(fi)),  for 

(5.6)  F(x,t)  =  bg(x)(Ug(x)-u(x,t) )  +  /p  h(xrs)ds  -  / g  B(x, t-s)u(x,s)ds 

e  1.2  2 

On  the  other  hand,  W  is  hounded  in  W  (0,T;  L  (J))).  In  order  to 

£ 

characterize  the  possible  weak  limits  W  of  W  (the  uniqueness  of  which 

will  follow,  as  usual  from  the  unique  solvability  of  the  limit  equation),  we 

shall  assume,  after  extracting  another  subsequence,  still  denoted  e,  that 
£ 

W  converges  weakly  to  some  W.  So,  going  to  the  limit  in  (5.3)  yields: 


The  main  task  is  to  find  the  relationship  between  W  and  u,  which  (5.4) 

should  yield.  Here,  because  the  problem  is  linear,  we  use  the  Laplace 

transform,  but  to  do  so  we  extend  the  problem  to  [0,+«)  in  time  as 

C  t  £  £ 

follows:  for  t  >  T,  extend  8  by  6  (x,t)  =  8  (x,T),  h  and  y  by  zero 
c  e 

(so  c  (x,t)  =  c  (x,T) ) ,  and  by  theorem  (4.20)  which  applies  to  any 

e 

interval  [0,T.j],  u  exists  for  all  t  >  0,  but  for  t  >  T,  the  problems 

become  simpler,  as  seen  from  (2.1): 

e 

bE  +  8e(T)uE  -  div  A(cftVue-Y*Vu£)  *  0  . 

u  at  x  u 

A  detailed  analysis  of  estimates  (4.6),  (4.13)  shows  that  in  this  particular 

case,  |uE(t)|  grows  at  most  exponentially  in  t  with  a  rate  uniform  in 

H  (fl) 

e.  Therefore,  all  the  Laplace  transforms  considered  here  will  be  convergent 
at  least  in  some  complex  right  half-plane  Re  X  >  A^. 

We  will  denote  by 

v( A )  ■  e  Atv(t)dt  for  Re  A  >  AQ 
and  (5.3),  (5.4)  yield  ( 5.8 )  ^  below  since  the  gradient  operator  in  x 
commutes  with  the  Laplace  transform. 

I  WE(x,A)  -  AE(x)  c£(x,A)  V(ue(x,A)) 

(5*8>e  { 

^  -  div  WE(x,A)  ■  FE(x,A) 

For  fixed  A,  (5.8  )£  is  just  the  homogenization  problem  for  an  elliptic 
operator  with  complex  coefficients. 

Upon  inspection  of  (5.5),  one  sees  that  for  every  t  >  0,  Fe  converges 

to  F  in  C([0,t]j  H  (Q))  but  that  F  grows  in  H~  at  the  same  rate  as 

0  A  0 
u  does.  Consequently,  for  each  A  with  Re  A  >  XQ,  F  (A)  converges  to 

F(A)  in  H  (0;  C).  Similarly  ug(A)  converges  to  u(A)  weakly  in 

H.l(fl;  C).  The  sesquilinear  form 


is  continuous  coercive  on  Hgtfli  C)  under  the  hypothesis 

A 

(5.10)  Re  c£(x,X)  >  Pq(X)  >  0,  which  we  will  check  later. 

A 

Indeed,  c  is  bounded  on  fl  for  each  X  such  that  Re  X  >  0  so  a 
E  £ 

is  continuous  and  for  such  X ' s , 

Re  u,u)  -  Jn  Re  ce(x,X )AC(x) | Vu(x) | 2dx  > 

A  M 

>  u  fa  Re  ce<x,X)|Vu|  dx 
e 

since  A  is  symmetric  real  coercive.  Incidentally,  another  proof  of 

e 

existence  for  the  solution  u  is  thus  obtained  by  applying  Lax-Mi 1 gram 1 s 
theorem  for  u  . 

One  can  now  apply  a  compactness  result  for  complex  homogenization  (see 
for  example  Sanchez-Palencia  [1],  Murat  [1]  or  Bensoussan-Lions- Papanicolaou 


[1]).  For  each  X  with  Re  X  >  XQ,  there  is  a  matrix- valued  function 

A  A 

D(x,X)  (independent  of  F  and  W)  such  that  (5.8>e  implies  at  the  limit 


€  ♦  0: 


(5.11) 


-  div  W  -  F 

W(x,X)  -  D(x,X )7u(x,X)  . 


In  order  to  apply  the  inverse  laplace  transform  to  (5.11),  all  that  is 
needed  is  that  D(x,X)  be  analytic  in  X  with  at  most  polynomial  growth  at 
|  X  f  ♦  •,  in  which  case  it  is  the  Laplace  transform  of  a  distribution  of 


finite  order  in  t,  denoted  E(x,t).  That  D  is  analytic  is  a  mere 


consequence  of  the  fact  that  it  is  a  limit  of  a  sequence  of  analytic  functions 
of  X,  the  limit  being  locally  uniform.  From  the  uniform  boundedness  of  cC 
and  AC,  one  can  conclude  that  D(x,X)  is  bounded  by  a  multiple  of 
(Re  X)  \  Consequently,  E(x,t)  is  a  bounded  distribution  of  order  not  more 


than  2,  on  [0,-h*[,  with  values  in  the  cone  of  bounded  measurable  symmetric 


square  matrices  on  ft. 

He  now  check  that  (5.10)  holds:  integration  by  parts  in 

(5.12)  Re  c{.(x/X)  ■  fQ  e  cos(Im  X  t)cg(t)dt 

gives 

Re  c£(x,X)  ■  JQ^^  -  cos(t  Im  X))e  **  Ncjl  -  2  Re  Xc^.  +  (Re  X)2c£)dt  . 
Since  (1  -  cos(t  Im  X))  and  -c^  are  nonnegative  functions  and  c”  is  a 
nonnegative  measure, 

(5.13)  Re  c£(x,X)  >  /^*°  (1  -  cos(t  Im  X))e  **  *(Re  X)2c£(t)dt  » 
combining  (5.12)  and  (5.13)  one  gets 

2 

.  _  .  . .  _  .  ,  .  .  ( Re  X )  -t  Re  X  .... 

(5.14)  Re  c  (x,X)  >  — 1 - L—r  J  e  c  (t)dt  / 

e  1+(Re  XT  0  e 

*  Pjg  \ 

But  c  (t)  >  o  implies  with  (5.14)  that  Re  c  (x,X)  >  a  - -  which 

e  6  1+(Re  \) 

implies  (5.10).  Finally,  we  have  proved  the  following  theorem: 

£  £  £  £  £ 

(5. 15)  Theorem:  Let  bQ,  c  ,  0  and  a  “  A  (linear  case)  satisfy 

hypotheses  (4.4)  with  a  and  k  independent  of  e.  There  exists  a  sequence 

e  converging  to  zero,  functions  bn(x),  $(x,t)  and  a  distribution  E(x,t) 
n  e 

such  that  the  solution  u  n  of  the  corresponding  problem  (4.1),  (4.2) 

2 

converges  in  C([0,t];  l  (ft))  to  the  solution  u  of 

(5.16)  bQu'  -  div((E*Vu)')  -  G(u)  . 

e  e 

bQ  and  0  are  the  weak  limits  of  bQ  and  0  ,  and  E  is  obtained  via 

its  Laplace  transform  D(x,X)  which  is  the  complex  elliptic  homogenization  of 

*  e_ 
c£  (x,X)A  (x).  ■ 

n 

Remark:  Even  in  the  case  of  periodic  problems,  where  there  are  explicit 
formulas  for  D  it  is  not  known  whether  t E(t,x)  is  in  some  appropriate 
sense,  a  convex  decreasing  function  of  t,  not  even  whether  it  is  a  function 

of  t,  as  one  would  suspect.  This  is  one  of  the  problem  left  open  in  the 
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theory,  the  other  one  being  the  homogenization  of  the  non-linear  case  (where 


the  Laplace- transform  cannot  be  used) . 
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